Colored Noise in Quantum Chaos 
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We derive a set of spectral statistics whose power spectrum is characterized, in the case of chaotic 
quantum systems, by colored noise 1/f , where the integer parameter 7 critically depends on the 
specific energy-level statistic considered. In the case of regular quantum systems these spectral 
statistics show l/f +1 noise. 
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Quantum chaos is the study of quantum systems which 
are classically chaotic [1]. There is a conjectured re- 
lationship, backed up by numerous examples, between 
the energy level fluctuation properties of a quantum sys- 
tems and the dynamical behavior of its classical analog. 
In particular, Berry and Tabor showed that classically 
integrable systems are characterized by Poissonian spec- 
tral fluctuations [2], while Bohigas, Giannoni and Schmit 
[3] suggested that classically chaotic systems can be de- 
scribed by Gaussian ensembles of random matrix theory 
(RMT). In the past years many energy-level statistics 
have been proposed to characterize the spectral fluctua- 
tions; among them there are the nearest-neighbour level- 
spacing distribution P(s), the spectral rigidity Aq(L) and 
the spectral form factor K(t). Recently, Relano et al. 
[4] have introduced the energy- level statistic S n , showing 
that chaotic quantum systems display 1/f noise in the 6 n 
statistic, whereas integrable ones exibit 1/f 2 noise [4,5]. 

In this Brief Report we introduce a new set of spectral 
statistics which contains the spectral statistic S n as a par- 
ticular case. We show that for chaotic quantum systems 
these statistics are characterized by colored noise 1/f, 
where the integer parameter 7 critically depends on the 
specific energy-level statistic considered. Moreover, we 
find that in the case of regular quantum systems these 
spectral statistics display l// 7+1 noise. 

Let us consider a quantum system with a discrete set 
of unfolded levels E n , such that the energy is measured in 
units of the mean level speacing. The staircase function 
N(E), which gives the number of energy levels up to the 
energy E, can be written as 



N(E) = N(E)+N osc (E) 



(1) 



where N(E) — E is the averaged number of levels and 
N osc (E) is the fluctuating part. N osc (E) is supposed to 
belong to a universality class, which should only depend 
on the integrability or chaoticity of the classical analog 
[1] . The Madrid group of Relano et al. [4] has analyzed 
these fluctuations by using the spectral statistic 8 n , de- 
fined as 



N osc (E n ) 



(2) 



In particular, the Madrid group has investigated the be- 
havior of the power spectrum S(k) of the discrete series 
5 n , given by 



S M (k) = |4 1 2 , 
where Sh is the Fourier transform of 8 n 



M 



5 n exp 



ikn 



M 



(3) 



(4) 



and M is the size of the series. By studying spectra of 
atomic nuclei at high energies and also spectra of Gaus- 
sian ensebles of RMT, the Madrid group has found the 
power law (£m(&)) ~ For Poisson spectra it has 

found instead (S^fTf)) ~ 1/fc 2 , as expected for indepen- 
dent random variables. Finally, the following conjecture 
has been suggested: the energy spectra of chaotic quan- 
tum systems are characterized by pink noise 1/f, in con- 
trast to the Brown noise 1/f 2 of regular systems [4]. 

The Madrid group has analytically proved its conjec- 
ture on the basis of the RMT [5] . Independently, Robnik 
has obtained the same result [6]: under the conditions 
M ^> 1 and k/M <C 1, the averaged power spectrum of 
the energy-level statistic S n is given by 



(S M (k)} 



for chaotic systems 



Try- for integrable systems 



(5) 



where [3 depends on the symmetry of the Gaussian en- 
semble: = 1 for the Gaussian Ortogonal Ensem- 
ble (GOE), = 2 for the Gaussian Unitary Ensemble 
(GUE), and j3 = 4 for the Gaussian Symplectic Ensemble 
(GSE). Thus, although the origin of 1/f noise in many 
complex systems is still an unsolved problem [7] , the ori- 
gin of 1/f noise in the spectral fluctuations of chaotic 
quantum systems has been easier to understand because 
of the mathematical tractability of RMT. It is important 
to stress that Eq. (5) is valid if the condition k/M <g; 1 is 
fullfillcd. In fact, as shown in [5] the Eq. (5) is the power- 
law approximation, obtained via Taylor expansion, of a 
more complex formula. In Fig. 1 we plot both the exact 
curve derived in [5] and the power-law curve of Eq. (5). 
Figure 1 shows that (iSjvf(fc)) presents deviations from 
the power-law behavior for k close to the Nyquist num- 
ber k c = M/2. These deviations are due to the finite 
size M of the series S n and the figure shows that they 
are negligible for k/M < 10~ 2 , where the relative error 



1 



is below 1%. In practice, because of the log-log scale, the 
differences are clearly visible only for k/M > ICP 1 . 




FIG. 1 (color online). Average power spectrum 
of the spectral statistic S„. Comparison between the exact 
formula of Ref. [5] (solid line) and its power-law approxima- 
tion (dotted line), given by Eq. (5). The dashed line is the 
relative error of the two curves. Note the different scale of 
the Poisson ensemble (integrable systems) with respect to the 
Gaussian ensembles (chaotic systems). 

The proof of the nice result of Eq. (5) is based on 
the relationship between Sjwik) and the widely studied 
spectral form factor of the density of levels [5,6]. Here, 
we apply the same methodology of Ref. [5,6] to introduce 
and analyze a new set of spectral statistics, which con- 
tains 5 n as a particular case. First we define the following 
set of functions 



osc [ ' dE<* ' 



(6) 



where the parameter a is a positive integer. This fam- 
ily of functions can be extended also to negative integer 
values of a by setting 



N^ al) (E) = [ dx x f 1 d 



dx\ a \N osc (x\ c 



(7) 

For a — one recovers the oscillating part N osc (E) of the 
staircase function, while for a = 1 one has the oscillating 
part Posc{E) of the density of levels. In analogy with 
Eq. (2) we now introduce the following set of spectral 
statistics: 



8^ = NtKE n ) 



(8) 



spanned by the integer number a. Obviously, for a = 
one has the spectral statistic 5 n of Eq. (2). 



By applying the formula of integration by parts, it is 
straightforward to prove that the Fourier transform of 
N osc (E), given by 



N osc (t) 



+ 00 



N oac (E)exp(-iEt) dE , 



(9) 



is related to the Fourier transform of Nosc(E) by the 
simple formula 



N osc (t) = N^ c {t)-^- 



(10) 



Under the conditions M ^> 1 and k/M <C 1 the discrete 
Fourier transform of S n can be obtained from the Fourier 
transform of N osc (E). It is given by 



5k = N osc(jj) 



(11) 



Taking into account Eqs. (3), (8), (10) and (11), the 
relationship between the power spectrum of 5 n and the 
power spectrum of 5 n °^ reads 



S M (k) 



(12) 



Finally, by using the Eq. (5) one finds that the average 

power spectrum of the spectral statistic 5^ satisfies the 
following formula 



2 Af 1 



M 2 ~ 2 ° 



for chaotic systems 



for integrable systems 



(13) 



This formula is the main result of the paper. For a = 
it reduces to Eq. (5). In general, for a > the aver- 
age power spectrum of 5^ is not divergent, it is instead 
divergent for a < 0. 

Among the new spectral statistics, for numerical pur- 
poses the most simple ones are 5n ^ and 5$. The statis- 
tic Sn is given by 

5 n ~ 1] = N osc (E) dE , (14) 

J — oo 

and its average power spectrum is 

for chaotic systems 

for integrable systems 



2 M 
,3~F~ 



(15) 



iU 1 
T 1 " 



Thus, S n shows black noise l// 3 for chaotic systems 
and black noise l// 4 for regular systems. The statistic 



sr(i) 



is instead given by 

_ dN osc (E n ) _ 
n ~ dE " 



:{E n ) 



(16) 



2 



Its average power spectrum is nothing else than the dis- 
crete spectral form factor of the density of levels, which 
shows blue noise / for chaotic systems and white noise 
/° = 1 for regular systems. 



satisfactory way on the basis of random matrix theory. 

The author thanks J.M.G. Gomez, A. Relano, E. 
Faleiro and M. Robnik for many useful discussions. 
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FIG. 2 (color online). Average power spectrum (S^'(fc)) 



of the spectral statistics dh ^ an d S^ 1 ' . Power laws for inte- 
grable systems (Poisson) and for chaotic systems with differ- 
ent classes of symmetry: GOE, GUE, and GSE. 

For the sake of completeness, in Fig. 2 we plot the 
power spectrum of spectral statistics dn (upper panel) 

and SrP (lower panel). Figure 2 shows that different 
power laws can be easily distinguished by choosing the 
appropriate scale interval in the log- log plot. The new 
spectral statistics we have introduced are a simple and 
useful tool to analyze the energy levels of quantum sys- 
tems and determine their integrability or chaoticity. 

An open problem is the behavior of (S^\k)} for sys- 
tems in the mixed regime between order and chaos. Re- 
cently we have analyzed the order to chaos transition 
in terms of the power spectrum (5m(&)) by using the 
Pascal's snail (Robnik billiard) [8]. We have numerically 
found a net power law l// 7 , with 1 < 7 < 2, at all the 
transition stages. We have suggested that the exponent 7 
is related to the chaotic component of the classical phase 
space of the quantum billiard, but a theoretical explana- 
tion of these numerical results is still lacking [9] . 

In conclusion, we have shown that colored noise l// 7 
characterizes the spectral fluctuations of quantum sys- 
tems. The presence of colored noise in fluctuating phys- 
ical variables is ubiquitous. This kind of noise has been 
detected in DNA sequences, quasar emission, river dis- 
charge, heartbeat, among many others. Despite this 
ubiquity, there is no universal explanation about this 
phenomenon, which is a generic manifestation of com- 
plex systems. Nevertheless, here we have shown that 
for integrable and chaotic quantum systems the colored 
noise of energy-level fluctuations can be explained in a 
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